One-dimensional semiconductor in a polar solvent: Solvation and low- frequency 

dynamics of an excess charge carrier 
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Due to solvation, excess charge carriers on Id semiconductor nanostructures immersed in polar 
solvents undergo self-localization into polaronic states. Using a simplified theoretical model for 
small-diameter structures, we study low-frequency dynamical properties of resulting Id adiabatic 
polarons. The combined microscopic dynamics of the electronic charge density and the solvent leads 
to macroscopic Langevin dynamics of a polaron and to the appearance of local dielectric relaxation 
modes. Polaron mobility is evaluated as a function of system parameters. Numerical estimates 
indicate that the solvated carriers can have mobilities orders of magnitude lower than the intrinsic 
values. 

PACS numbers: 71.38.-k, 73.63.-b, 31.70.Dk 



INTRODUCTION 

One-dimensional (Id) semiconductor (SC) nanostruc- 
tures like nanowires and nanotubes are fascinating ob- 
jects offering a promise for various applications, includ- 
ing electronics, energy harvesting and sensors. The 
properties of excess charge carriers are fundamental for 
these applications. A customary picture used to describe 
charge carriers in nanowires 'll| and nanotubes '2!| is that 
of band electron states with very high, nearly ballistic [3| , 
intrinsic mobilities. Many applications, however, such as 
photoelectrochemical energy conversion 4, S' , feature SC 
nanostructures in contact with polar liquid environments. 
We have recently pointed out [g, |7| that excess carriers 
on Id SCs immersed in the sluggish polar medium can 
drastically change their nature. Due to the long-range 
Coulomb interaction with the slow (orientational) po- 
larization of the medium, a charge carrier (an electron 
or a hole) gets solvated forming a localized electronic 
state surrounded by a self-consistent dielectric polariza- 
tion pattern as shown schematically in Fig. [TJ In phys- 
ical parlance, self-trapped states formed by this mecha- 
nism are traditionally called polarons, extensively stud- 
ied in 3d ionic crystals and polar SCs [8|, |9[. Another 
related 3d notion is that of the solvated electron in liq- 
uids [ll[ll[ll. 

The formation of Id polarons is expected to have an 
impact on various processes such as electron-hole sepa- 
ration from bound excitonic states [6|j7| , charge transfer 
reactions, local optical absorption [l3J and charge car- 
rier transport. As we discuss the low- frequency polaron 
dynamics in this paper within a simplified theoretical 
framework, it will be particularly illustrated that the mo- 
bilities of the solvated carriers can be orders of magnitude 
lower than the intrinsic values. We note that the system 
we consider appears an interesting model realization in 
the general context of quantum particles interacting with 
a dissipative environment [ij] . While we restrict our dis- 
cussion to rigid SC nanostructures with good intrinsic 
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FIG. 1: Schematically: Solvated excess electron charge den- 
sity (thick line) along the Id semiconductor nanostructure 
which is surrounded by a pattern of the orientational polar- 
ization of solvent molecules. 



electronic properties and favoring the formation of con- 
tinuum adaibatic polarons, the generic aspects apply to 
other systems as well. Noteworthy, Basko and Conwell 
[15J emphasized a conceptually similar formation of hole 
polarons in DNA in aqueous electrolyte solutions and dis- 
cussed the effect of water drag on their mobility [16[ . 

The parametric scales for the spatial extent and the 
binding energy of the 3d polaron are established by the 
3d Coulomb center problem where, for a carrier of charge 
q and effective mass m in the uniform medium of effective 
dielectric constant e*, the corresponding 3d Bohr radius 
and Rydberg are [30| 



as 



e*h jvaq^ , Ry = q^ /2e*aB- 



(1) 



As a result of the confinement of charge carrier motion to 
the Id SC, the effects of Coulomb interactions are gener- 
ally known 17[ to get amplified, and, in fact, the binding 
energy E^° of the confined polaron can increase substan- 
tially in comparison with the 3d binding [6|,|7|. Our inter- 
est is in the strong confinement regime: as ^ R {R is the 
transverse radius of the Id structure), when the resulting 



longitudinal spatial extent of the polaron is much larger 
than R (see Fig. [2jb)) and the self-localized charge car- 
rier can be adequately described by the Id wave function 
ip{x), X being the coordinate along the structure axis. 

We consider the polaron to be formed due to the in- 
teraction of an excess carrier with the slow polarization 
component of the Debye solvent described by the fre- 
quency w-dependent dielectric function 



e{u;) 



1 — iujTo ' 



(2) 



for which the effective dielectric constant in Eq. ^ would 



be determined from the well-known ^ 
ship 

1/e* = 1/eoo - 1/e., 



Ill, [IJI relation- 



(3) 



with es ^ eoo for typical solvent values of static ts 
and high-frequency too constants [18| . The characteristic 
time Tu in Eq. ^ is the Debye (transverse) relaxation 
time, for the problem at hand, however, the solvent re- 
sponse occurs at generally much shorter time scales cor- 
responding to the longitudinal relaxation time 12l.ll8Lll9l| 



TL 



■ro- 



(4) 



Typical solvents are charaterized by a wide range of r^ 
ranging from fractions to tens of ps [y, ll8[. With the 
estimated polaron binding energies E\° on the order of 
0.1 eV, the dynamic adiabaticity is ensured by Ti/r/, ^ 
-E^° . For room temperatures T, on the other hand, one 
can also have fc^T ^ _E^° satisfied, thus making thermal 
population of higher energy electronic states negligible. 

The aim of this paper is to elucidate the response of 
such Id polarons to weak perturbations of low- frequencies 
w: fiLo ^ E\° , which do not cause transitions between 
different electronic levels. This, in particular, encom- 
passes such basic phenomena as polaron drift caused by 
an applied field and polaron diffusion caused by the di- 
electric fluctuations of the surrounding solvent. We stress 
that our consideration is limited to a simplified theoreti- 
cal model whose main purpose is to address elements of 
the essential physics. We will show how the underlying 
combined microscopic dynamics of the electronic charge 
density on the Id SC and the solvent results in point- 
particle- like Langevin dynamics of the polaron center. 
The derived mobility and diffusion coefhcient, as should 
be expected 1J|, satisfy the Einstein relation. In ad- 
dition, we demonstrate the existence of local relaxation 
modes around the polaron whose relaxation times are 
somewhat longer than bare solvent's r^. These modes 
involve polaron shape variations and presumably could 
manifest themselves in responses to time-dependent fields 
as well as in collisions. 



GENERAL ADIABATIC POLARON 
RELATIONSHIPS 

In the adiabatic picture, the excess charge carrier 
instantaneously responds to the electrostatic potential 
0(a;, i) that is slowly changing in time by staying in the 
lowest eigenstate of the stationary Schrodinger equation: 

Hi:{x,t)^E^{x,t), H = - — —^ + q(t)ix,t), (5) 

Zm ax'^ 

where time t thus enters only as a parameter. The cor- 
responding time-dependent Id charge density is 



p{x,t) = q\ipix,t)\'' 



(6) 



where q is the charge of the carrier. 

The electrostatic potential (j)[x,t), in turn, can be rep- 
resented as a sum of two contributions: 



0(x,t) = (j)i{x,t) +(f)six,t), 



(7) 



where the first term represents the potential "indepen- 
dent" of the charge carrier itself, that is, the potential 
that would be there if the charge carrier was absent - as 
exemplified by applied fields and dielectric fluctuations 
of the environment. The second term in Q is the part 
that is induced by the carrier charge density ([6]) in the 
past via the time-dependent medium response: 



dt' / dx'Q{x' ,t')p{x — x' ,t ~- 1'). (8) 



<t>six,t) 



Equation ([8]) assumes the translationally invariant re- 
sponse so that 

(j)sik,Lu) = Q{k,uj)p{k,ijj). 

in Fourier space. The response function Q is determined 
from the solutions of the corresponding 3d electrostatic 
potential problem using the appropriate geometry of el- 
ementary charge distributions and of the dielectric en- 
vironment. Let us denote such a solution derived for a 
particular value c(u}) of the medium dielectric function 
as V{k,e{uj)). One, of course, recognizes that Q should 
not include the instantaneous Coulomb interaction of the 
charge carrier with itself; in the context of the Debye sol- 
vent ([2|) , the medium for this case is characterized by Coo ■ 
In these terms. 



Q{k,uo)^V{k,e{uj)) ~V{k,e^), 



(9) 



the relevant response is thus due to the slow (orienta- 
tional) component of the dielectric polarization of the 
medium. Explicit examples of such calculations for the 
nanotube geometry with the account of tube's own po- 
larizability can be found in Ref. [6|. When the dielectric 
polarization effects are well captured by the model of a 
uniform 3d medium, which is a good approximation in 



the strong confinement regime [a, |20|, the response Q 
is separable into a product of the spatial and temporal 
components (all bare Coulomb interactions are simply 
decreased by a factor of the inverse dielectric constant). 
We will write this as 



g{x,t)^-A{t)G{x) 



g{kuj)^-a{uj)g{k). (10) 



To gain more physical insight and obtain explicit simpler 
relationships, we will restrict our attention in this paper 
to this regime. In accordance with Eq. ([9|), the temporal 
part for the Debye solvent ([2|) can be described as 



a{uj) = 



1 



1 



"0 



1 



e{Lu) 1 — iuJTL ' 



ao 



(11) 



The relaxation time T]^ that determines the pole of the 
inverse dielectric function as well as the familiar com- 
bination® naturally arise here. As is well known 



I2I llSl . Iig | , td characterizes the dielectric relaxation in 
circumstances controlled by electrostatic potentials (elec- 
tric fields E) while tl is appropriate in circumstances 
controlled by charge densities (electric inductions D). 
Evidently, the latter is the relevant situation here, and we 
will use T = Ti to simplify the notation in what follows. 
In the time domain, the response function is therefore 



, / X f duj 
^^'^ = J ^ 



— *c,(^)^^e^*/^« 



e(t), (12) 



where is a step function. 

The spatial part G{x) of the response p^ is deter- 
mined by the specifics of the transverse charge distribu- 
tion in our Id SC. For the illustrative calculations below, 
we will use the case of the nanotube geometry where (see, 
e.g., Refs. [i,|2l|) 



m = 2Io{kR)Ko{kR), 



(13) 



with Iq and Kq being the modified Bessel functions and 
R the tube radius. Equation (|13p essentially corresponds 
to the Fourier transform of the Coulomb interaction be- 
tween two charged rings of radius R. 

Taking advantage of the particular simple form P^ . 
the temporal evolution of the potential (j){x,t), as per 
Eqs. dH) and ([8|), can be conveniently converted into the 
differential (over time) equation: 

d(t>{x,t) 4>{x,t) d4>i{x,t) (i)i{x,t) 



dt 



dt 



T 



dx'G{x~x')p{x',t). (14) 



THE STATIONARY SOLUTION AND SMALL 
PERTURBATIONS AROUND IT 

One first discusses the dynamics in the absence of 
the external potential (j)i{x,t), and the starting point, of 



course, is the consideration of the ground-state stationary 
polaron problem. In this basic case, the resulting charge 
density and potential are time- independent, we denote 
their spatial distributions as p{x) and (^(x), respectively. 
Their relationship is determined by the static medium 
response, or a(0) = ao in Eqs. pI7)l . (fTT|) . That is also 
immediately seen from Eq. (fT4|) : 



<j){x — xo) — — ao / dx'G{x — x') p{x' — xq), (15) 

where we explicitly emphasized the translational invari- 
ance by indicating that the stationary solutions can be 
centered around arbitrary points xq. (The negative sign 
in the r. h. s. of Eq. (fT5|) and in other expressions below 
signifies that a positive charge density leads to a sta- 
bilizing negative potential, and vice versa.) As p{x) is 
determined ([5]) by the ground-state charge carrier wave 
function ipo{x) obeying Eq. ([5]), the latter wave func- 
tion is then found from the self-consistent non-linear 
Schrodinger equation 

~£^^ ~ ^ Jd^'Gix-x')mx')\'^{x) ^ E^{x) 

as corresponding to its lowest eigenvalue E — Eq. This, 
of course, is the same equation one would derive from the 
optimization of the adiabatic system energy; it results 
in self-localized (polaronic) solutions as we discussed in 
Refs. |6|,l7|. 

With the ground-state polaron solution in place, one 
can now consider the effects of small potential perturba- 
tions (j){x^t): (j){x^t) = (l){x) + 4>{x,t), in the stationary 
Schrodinger problem ([5]): H = H + V, 



H 



2m dx^ 



+ q(f){x), V^q4>{x,t). 



The overall approach here is very similar to used in previ- 
ous studies of small perturbations around 3d [23| and Id 
23I ]2M, ]2M optical polarons. The standard lowest-order 
perturbation theory with respect to V then results in the 
perturbed electron density p{x, t) — p{x) + p{x, t) with 



p{x,t) = — I dx'L{x,x') (j){x' ,t), 



(16) 



where the Hermitian kernel 

tl;^{x)tl;^{x)il;*{x')'ipo{x') + c.c. 

i>0 



Lix,x')=q'J2- 



Ei — Eq 

is built out of the solutions of the unperturbed problem 

H^^{x) = Eiipi{x). 



The density perturbation (fT6|) should self-consistently de- 
termine the potential perturbations ([8]) at future times, 
that is: 

4>{x,t)= I dt' I dx'g{x',t')p{x-x',t-t'). (17) 



Together, Eqs. (fTBl) and ((T7|) determine the eigenmodes 
of the un-driven dynamics around the stationary polaron 
in the case of a general response G- In the separable 
case (|10p , the problem of finding eigenmodes evidently re- 
duces to a simpler generalized eigenvalue problem which 
we write in the form of 



Pn{x) = - dx'L{x,x')(j>n{x'), (18a) 

(f>n{x) = — A„ao / dx'G{x — x') pn{x'). (18b) 



In practical terms, the generalized problem (|T8|) is easily 
converted to a normal eigenvalue problem and readily 
solved in Fourier (wave vector) space. The eigenf unctions 
of this problem form a complete orthogonal set. To save 
space on notation, here we choose to normalize this set 
without extra coefficients as 



/ 



dx Pmix) (pnix) = S„ 



(19) 



(all functions can be chosen real) . The eigenvalues A„ > 1 
(n > 0), in turn, determine the eigenfrequencies w„ from 



a{u!n) = A„ao. 



(20) 



Unhke the oscillating solutions in Refs. [H, [H, [H, [111 , 

for the Debye solvent with its dissipative dynamics pTj) . 
Eq. ([20|) leads to a set of relaxation times r„ : 



UJ„ 



I 



1 - 1/A„ 



(21) 



The temporal evolution of the eigenfunctions of the dy- 
namical problem would therefore be "ordinarily" decay- 



ing: 



p„{x,t) oc Pnix)e */'^" 



Figure [^a) illustrates how these relaxation times t„ con- 
verge to the "bare" value of r upon the increase of the 
mode index n - for this illustrative calculation we used 
the tubular case (J13p with the confinement parameter 
Ub/R = 10. One immediately notices, however, that the 
dynamics has a very different character for the mode with 
Ao = 1 corresponding to o^o = 0. This is the so-called 
zero-frequency translational mode whose existence is a 
general feature of quantum-field models with non-trivial 
spatial solutions and is a consequence of the translational 
invariance [20|. As discussed above, the static polaronic 
solutions are degenerate with respect to the position of 
the polaron center. An infinitesimal displacement along 
X corresponds to the spatial derivative of the static solu- 
tions over X (denoted below with symbol '); one indeed 
easily verifies that 



p^{x)^^-'/^-p'{x), c^o{x)=l-^'^'{x) 
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FIG. 2: Illustration of local relaxation modes (tubular geome- 
try with ub/R = 10): (a) Relaxation times for the first n — 
to 8 eigenmodes; (b) Shapes of the spatial density variations 
pn{x) for the first 3 modes as indicated by values of the index 
n next to thin lines. These lines are arbitrarily scaled here to 
be better displayed. The thick line corresponds to the scaled 
static polaron density aBp{x)/q. 



satisfy Eqs. (fT5|) with Ao = 1. 
factor involves 



Here the normalization 



7= dx\p'{x)^'{x)\. 



(22) 



(23) 



As an example, Fig. [DJb) shows spatial patterns corre- 
sponding to the first 3 eigenmodes in comparison with 
the static polaron shape in terms of the charge density. 
Of course, all density variations satisfy J dx pn{x) = 0. 
Given the numerical nature of the solution, we have con- 
firmed that the translational mode {n — 0) precisely fol- 
lows the spatial derivative of the static polaron shape. 
The next mode (n — 1) can be qualitatively thought of 
as corresponding to the variation of the polaron width. 
The convergence pattern in panel (a) and the spatial vari- 
ations displayed in panel (b) of Fig. [2] are reminiscent of 
those discussed in connection with Id optical polarons 



As is well known [2a, [27 1, in the expansions of the 
arbitrary spatial patterns over the complete set of eigen- 
modes, the zero-frequency mode should be excluded in 
favor of the collective coordinate ^{t) of the polaron cen- 
ter position (polaron centroid). Such expansions for the 
potential and density variations then take the form of 

cl>{x,t) = 4>{x - ^{t)) + Y,an{t)Mx - m), (24a) 

n>l 

p{x,t) = p{x-^{t)) + J2<^n{t)pn{x-m), (24b) 
n>l 

where the the first terms correspond to the shapes of 
the static solutions that are centered around ^{t). The 
set of time-dependent coefhcients a„(t) (n > 1) and 
the centroid coordinate ^(i) then constitute complete co- 
ordinates of the system and can be conveniently used to 
describe its dynamics. 

Using Eqs. (|24|) in Eq. (fT4|) for the potential evolution, 
one converts the latter into 



and 



e 



'' + I]a„0'„ 



n>l 



E 

n>l 



(/)„ = F{x,t), 
(25) 



where the external perturbation part is 

d(j)t{x,t) (t)i{x,t) 



F{x,t) = 



dt 



(26) 



As in Eqs. ([M]) . all potential expansion functions in 
Eq. pS)) have the same argument {x — ^{t)) . 



EQUATIONS OF MOTION AND CENTROID 
DYNAMICS 



One derives the equations of motion for the set of co- 
ordinates {^(i),a„(i)} by projecting Eq. (^5)1 with each 
of Pn{x — £,{t)) (n > 0) and utilizing the orthonormality 
relationships p^ . Now note that we are interested in 
and have been developing our framework to study the 
lowest-order effects of small perturbations. Correspond- 
ingly, the system responses to such perturbations, the 
centroid velocity ^{t) and shape deformation coefficients 
a„(i), are small. Therefore, the second-order term in 
Eq. ([^S)) (containing products of ^ and a„) can safely 
be dropped. The result is then a set of fully decoupled 
equations of motion: 

7^ = /, (27) 

an+an/Tn = fn, n > 1, (28) 



where 



/(e,t)= dxp'ix)F{x + ^,t) (29) 



/«(C,t) 



dxpn{x)F{x + (,t). 



Equations ([77|) and (pS)) evidently allow one to easily 
study various dynamic responses. Here we will concen- 
trate on the centroid dynamics. 

The external (perturbation) potential <j)i[x,t) can be 
represented as a sum of two contributions: 

(j)i{x,t) = (j>a{x,t) +(j)f{x,t), 

where 4>a{x,t) is due to the applied field while 4)f{x,t) is 
due to the "bare" fluctuations of the dielectric polariza- 
tion (here, of course, only its slow, orientational, compo- 
nent is relevant). We will correspondingly use subscript 
indices a and /. 

Let us consider first the motion of the polaron under 
the action of a constant applied potential 4'a{x) that can 
slowly (on the scale of the polaron size) change in space. 
From Eqs. (pS)) and (P^ it immediately follows that 



faiO = / dxp'ix) 



Mx + 



dxp{x) 



= dxp[x) — — 

9^(0 



•a{x+£,) 



di 



(30) 



where 8 is the electric field acting on the polaron. 

On the other hand, the fluctuating contribution //(^, t) 
is a random process whose correlations have to be de- 
termined from the correlations of (j)f{x,t) subsequently 
using definitions (P^ and ([^5]) . The bare classical 
(fc^T ;» Ti/t) fluctuations are easily found by using the 
fluctuation-dissipation theorem for the specified Id ge- 
ometry, the result is 

{(t)f{x,t)4)f{x',t')) ^aokBTG{x-x')e-^'-*'^/\ (31) 

With Eq. (|26p . this translates exactly into 

2aokBT 



■G{x-~x')5{t~t'). (32) 



{Ff{x,t)Ff{x',t')) 



Given the time-uncorrelated character of Eq. ([32|) . the 
relevant correlations of f f are those for the same centroid 
coordinate: 

{fM,t)ff{U))^'^-:^^^5{t-t'). (33) 

T 

To derive Eq. ((331) from Eqs. ^^ and ([32]), we have also 
used Eq. ([T5|) in the following: 



dxdx'p'{x)p\x')G{x-x') 



dx 



P'{x)cj)'{x) 7 



Oto 



ao 



Introducing the random velocities ri{t) = ff{t)/^, the 
equation of motion (|27p for the polaron centroid under 



the action of a constant field can be conveniently rewrit- 
ten as 



i^ fiqS + ij{t), 



where the mobility 



/i = 1/7' 



(34) 



(35) 



is defined with respect to the dynamical force q£ while 
the random process ri{t), as per Eq. (|33p . satisfies 



with the diffusion coefficient 

D = kBT/-fT. 



(36) 



(37) 



Equation (I34p is the standard overdamped Langevin 
equation [12] for a point particle as consistent with our 
assumption of the field having practically no spatial vari- 
ation on the scale of the polaron size. We emphasize that 
ingredients of Eq. p4p have been derived from the un- 
derlying microscopic dynamics. Our mobility (j35p and 
diffusion coefficient ([57]) manifestly obey the Einstein re- 
lation 

D = fxksT, 

thus ensuring that the macroscopic dynamics of the po- 
laron in the steady state leads to the Boltzmann distri- 
bution. We note that expression ((35)) is similar in form 
to the one obtained from energy conservation considera- 
tions by Basko and Conwell 15|, [iGj but with the impor- 
tant quantitative distinction that our r = tl is the lon- 
gitudinal relaxation time rather than td they used (the 
difference that can easily exceed an order of magnitude). 

Our discussion of the constant applied field is 
straightforwardly generalized to time-dependent fields. 
As Eq. (PSJ) suggests, for fields varying as £{t) — 
£oexp{—iujt), the first term in Eq. ([M|) would get modi- 
fied, with 1/r being replaced by —{0;+ 1/r. In addition, 
time-dependent fields would be exciting polaron shape 
deformations ([28]) as we will discuss elsewhere. 

The calculation of the mobility ([35]) . within the stud- 
ied model, is thus reduced to the evaluation of Eq. ([23]) 
based on the ground-state of the polaron. Figure [3] dis- 
plays results of such a calculation for the tubular inter- 
action ([T3|) in the scaled form that separates the specific 
dependence of the mobility on the confinement parame- 
ter ubI R- Within the shown region of its variation, the 
behavior exhibited in Fig. [3] may be approximated by a 
power-law: 



(m/Mo) oc {as/R) 



-1.4 



(38) 



For the overall dependence of the polaron mobility on 
system parameters, one should recall that the reference 
value of 







a^l R 



FIG. 3: The polaron mobility in the constant applied field 
as a function of the confinement parameter as/R, here the 
mobility is measured in units of fio = a^/r^ Ry which itself 
can strongly depend on system parameters. 



is itself strongly parameter-dependent as per Eqs. ([T]) and 

O- 

Using Eqs. ([55)1 and ([55)) . one can establish how the 
mobility scales when some of system parameters vary. 
The effect of the radius R on the mobility is already ex- 
plicit in Eq. (|38p . For the dependence on the effective 
mass m, one finds an approximate scaling law 



11 (X m 



(40) 



The dielectric properties of the solvent are also impor- 
tant here: if one were to separate the effect of r^ on /i 
the remaining factor for the mobility would feature the 
dependence on e* as 



TL ■ ^J,cc (e*) 



2.6 



(41) 



Mo 



j/TiRy 



(39) 



With typical Cj ^ Eoo, e* — £00- To appreciate the ef- 
fect on the mobility itself, however, one is reminded that 
the longitudinal time tl (H]) is also proportional to Coo- 
Equations (PD|) and (|^T|) illustrate a quite understandable 
general trend: the more strongly bound the polaron is, 
then the lower its mobility would be. 



DISCUSSION 

One-dimensional semiconductor nanostructures in im- 
mediate contact with polar liquids is an interesting class 
of systems of relevance to many applications, particu- 
larly those involving redox reactions [J, |28| . Using sim- 
plified theoretical models, we have demonstrated [6J, |7| 
that the solvation of excess charge carriers due to the 
solvent polarization can cause carriers on Id SCs to self- 
localize. A combination of several factors, small intrinsic 
effective masses of the carriers, enhanced Coulomb ef- 
fects in Id (small-diameter structures) and the sluggish 



nature of the orientational solvent polarization, creates 
conditions in these systems that may be favorable for the 
resulting formation of Id continuum adiabatic polarons. 
The physical properties of the polarons are quite different 
from those of the band states. The energetic significance 
of the polaronic effect is emphasized by making compar- 
isons with much more studied Wannier-Mott excitons in 
Id SCs. When compared theoretically on an equal foot- 
ing, we have found [y, [Zl] that the polaron binding energy, 
E^° , can be a substantial fraction, roughly one-third for 
charge-conjugation-symmetric SCc, of the exciton bind- 
ing energy, E^^'^ , in the media with comparable €00 ■ 

Let us take a widely known example of single- wall car- 
bon nanotubes (SWCNTs). The importance of the exci- 
tonic effects in the optics of semiconducting SWCNTs is 
well established now, with binding energies E^^"^ exper- 
imentally measured in some tubes to be in the range of 
0.4-0.6 eV [2i,|30,|3H. That would suggest that polaron 
binding energies can reach values of 0. 1 eV and more, as it 
also follows from numerical estimates, reaffirming the po- 
tential importance of polarons. It would be pertinent to 
note that redox chemistry of CNTs has been deemed an 
"emerging field of nanoscience" [S^l and solvatochromic 
effects in CNTs are being intensely researched [33[; in- 
terestingly, even the first mapping of luminescence ver- 
sus absorption spectra of individual SWCNTs has been 
achieved in aqueous suspensions. |34l35| 



In an effort to establish experimentally testable signa- 
tures of polaron formation, in this paper we have stud- 
ied the low-frequency dynamics of Id large adiabatic po- 
larons. We have shown how the combined microscopic 
dynamics of the electronic charge density on the Id SC 
and the solvent results in macroscopic Langevin dynam- 
ics of the polaron centroid. The mobility ((35)) and diffu- 
sion coefficient ([37)) are thus derivable from the stationary 
polaron parameters and the solvent dynamic relaxation. 
We also demonstrated that new local dielectric relaxation 
modes develop in the presence of the polaron. These 
modes involve polaron shape deformations and may play 
a role in time-dependent responses and in polaron colli- 
sion events. 

It is instructive to make some numerical estimates 
based on the obtained results. If, for instance, one 
chooses a set of reasonable representative numerical pa- 
rameters: m = 0.05me (we being the free electron mass), 
e* = 3 and tl = I ps, then, by Eq. ((T|), the resulting 
as ^ 32 A, Ry ~ 76 meV, and, by Eq. ([39]), g ■ ^o - 1-3 
cm^/V-s. Using the data of Fig. [3] for qb/R = 6 as 
an example of a small-diameter structure, this value of 
liQ would then translate into 9 • /x ~ 0.3 cm^/V-s. At 
T=300 K this yields the value of the diffusion coefficient 
D ~ 8 X 10~^ cm^/s. Of course, with a wide range of 
variability of the relevant system parameters (suffice to 
mention values of t^), these tentatively estimated values 
can change substantially for various systems. Still, these 
numbers strongly indicate that the mobility of the sol- 



vated carriers can be orders of magnitudes smaller than 
the intrinsic carrier mobilities of lO'^-lO^ cm^/V-s dis- 
cussed for SC wires and tubes [l|, l2| . 

We conclude by reiterating that the illustrative calcu- 
lations presented in this paper are based on simplified 
generic models whose main role is to elucidate the essen- 
tial physics. While we do not expect substantial changes 
in the results valid for small-diameter nanostructures 
(the strong-confinement regime, when most dielectric 
effects are due to the surrounding medium [O, [Zl, 120]), 
more accurate calculations for specific systems, including 
a finite size of dielectric cavities, SCs band structure and 
polarizability, may be required for detailed comparison 
with prospective experimental data. 

We are grateful to V. M. Agranovich and E. M. Con- 
well for useful comments. This study was supported by 
the Collaborative U. T. Dallas - SPRING Research and 
Nanotechnology Transfer Program. 
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